Theory of the oscillatory photoconductivity of a 2D electron gas 
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We develop a theory of magnetooscillations in the photoconductivity of a two-dimensional electron 
gas observed in recent experiments. The effect is governed by a change of the electron distribution 
function induced by the microwave radiation. We analyze a nonlinearity with respect to both the dc 
field and the microwave power, as well as the temperature dependence determined by the inelastic 
relaxation rate. 
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Recent experiments have discovered [1] that the re- 
sistivity of a high-mobility two-dimensional electron gas 
(2DEG) in GaAs/AlGaAs heterostructures subjected to 
microwave radiation of frequency to exhibits magne- 
tooscillations governed by the ratio uj/uj Ci where oj c is the 
cyclotron frequency. Subsequent work [2-6] has shown 
that for samples with a very high mobility and for high 
radiation power the minima of the oscillations evolve into 
zero-resistance states (ZRS). 

These spectacular observations have attracted much 
theoretical interest. As was shown in Ref. [7], the ZRS 
can be understood as a direct consequence of the os- 
cillatory photoconductivity (OPC), provided that the 
latter may become negative. A negative value of the 
OPC signifies an instability leading to the formation of 
spontaneous-current domains showing zero value of the 
observable resistance. Therefore, the identification of the 
microscopic mechanism of the OPC appears to be the key 
question in the interpretation of the data [1-6] . 

A mechanism of the OPC proposed in Ref. [8] is 
based on the effect of microwave radiation on electron 
scattering by impurities in a strong magnetic field (see 
also Ref. [9] for an earlier theory and Ref. [10] for a 
systematic theory). An alternative mechanism of the 
OPC was recently proposed in Ref. [11]. In contrast to 
Refs. [8-10], this mechanism is governed by a radiation- 
induced change of the electron distribution function. Be- 
cause of the oscillations of the density of states (DOS), 
v(e), related to the Landau quantization, the correction 
to the distribution function acquires an oscillatory struc- 
ture as well. This generates a contribution to the dc 
conductivity which oscillates with varying lo/uj c . A dis- 
tinctive feature of the contribution of Ref. [11] is that it is 
proportional to the inelastic relaxation time T in . A com- 
parison of the results of Refs. [10] and [11] shows that the 
latter contribution dominates if r m 3> r q (where T q is the 
quantum, or single-particle, relaxation time due to impu- 
rity scattering), which is the case for the experimentally 
relevant temperatures. 

The consideration of Ref. [11] is restricted to the regime 
which is linear in both the ac power and the dc electric 



field. The purpose of this paper is to develop a com- 
plete theory of the OPC governed by this mechanism, in- 
cluding nonlinear effects. We will demonstrate that the 
conductivity at a minimum becomes negative for a large 
microwave power and that a positive sign is restored for 
a strong dc bias, as it was assumed in Ref. [7]. 

We consider a 2DEG (mass m, density n e , Fermi 
velocity vf) subjected to a transverse magnetic field 
B = (mc/e) lo c . We assume that the field is classically 
strong, w c T tr 3> 1, where r tr is the transport relaxation 
time at B = 0. The photoconductivity cr p h determines 
the longitudinal current flowing in response to a dc elec- 
tric field £dc, 3 ■ £dc = CTph^dc in the presence of a mi- 
crowave electric field £ u cos cot. The more frequently mea- 
sured [1-3,5,6] longitudinal resistivity, p p h, is given by 
p p h — p^yCph, where p xy ~ eB/n e c is the Hall resistivity, 
affected only weakly by the radiation. 

We start with the formula for the dc conductivity: 



(1) 



(T p h = 2 / efecr dc (e) [-d e f(e)} 



where /(e) is the electron distribution function, and 
o~dc{£) determines the contribution of electrons with en- 
ergy e to the dissipative transport. In the leading approx- 
imation [10,11], (Tdc(£) = °dc^ 2 ( e )' wnere ^( £ ) = v {z)l v Q 
is the dimensionless DOS, isq — m/2ir is the DOS per 
spin at zero B (we use ft — 1), and cr^ c = e 2 v a v\ / '2w 2 T tr 
is the dc Drude conductivity per spin. All interesting ef- 
fects are due to a non-trivial energy dependence of the 
non-equlibirum distribution function /(e). The latter is 
found as a solution of the stationary kinetic equation 
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where the ac Drude conductivity per spin is given by (we 
assume \u> ± Wc| r tr ^> 1) 
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On the right-hand side of Eq. (2), inelastic processes are 
included in the relaxation time approximation (more de- 
tailed discussion of the relaxation time T; n is relegated to 
the end of the paper), and /r(e) is the Fermi distribution. 
The left-hand side is due to the electron collisions with 
impurities in the presence of the external electric fields. 
The first term describes the absorption and emission of 
microwave quanta; the rate of these transitions was cal- 
culated in Ref. [11]. This term can be also extracted 
from the kinetic equation of Ref. [10]. The second term 
describes the effect of the dc field and can be obtained 
from the first one by taking the limit u> — > 0. 

Equation (2) suggests convenient dimensionless units 
for the strength of the ac and dc fields: 
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Note that V u and Qdc are proportional to T m and are 
infinite in the absence of inelastic relaxation processes. 

We consider first the case of overlapping Landau levels 
(LLs), with the DOS given by v = 1 - 2<5cos ^p, where 
8 = exp(— ir/uj c T q ) <C 1. Here r q is the zevo-B single- 
particle relaxation time, which is much shorter than the 
transport time in high-mobility structures, r q <C r tr (be- 
cause of the smooth character of a random potential of 
remote donors). The existence of a small parameter 8 
simplifies solution of the kinetic equation (2). To first 
order in 8, we look for a solution in the form 



f = f0+.fosc+O(5 2 ), f 
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We assume that the electric fields are not too strong 
[VUu/Tf <C 1 and Q dc (uj c /T) 2 <C 1], so that the 
smooth part /o(e) is close to the Fermi distribution /t(e) 
at a bath temperature T 3> ui c ; otherwise, the tempera- 
ture of the electron gas is further increased due to heat- 
ing. Smooth functions /0.1(e) change on a scale of the 
order of temperature. We obtain 
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and substitute Eq. (6) into Eq. (1). Performing the 
energy integration in Eq. (1), we assume (in confor- 
mity with the experiment) that T is much larger than 
the Dingle temperature, T 3> 1/2ttt q . The terms 
of order 8 in Eq. (1) arc exponentially suppressed 
8 J ded £ f T cos oc <5 exp(-27r 2 T/u c ) <C S 2 and can 
be neglected. The leading uj dependent contribution to 
cr pn comes from the S 2 term generated by the product of 
<9 e /osc(£) °c $ cos ^zr an d the oscillatory part —28 cos ^p 
of v{e). This term does survive the energy averaging, 
-/ ded e f T cos 2 2p ~ 1/2. We thus find 
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Equation (7) is our central result. It describes the pho- 
toconductivity in the regime of overlapping LLs, includ- 
ing all non-linear (in £ w and £dc) effects. Let us analyze 
it in more detail. In the linear-response regime (£dc — ► 0) 
and for a not too strong microwave field, Eq. (7) yields a 



correction to the dark dc conductivity a dc 
which is linear in the microwave power: 
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in agreement with Ref. [11]. It is enlightening to compare 
Eq. (8) with the contribution of the effect of the ac field 
on the impurity scattering [8-10]. The analytic result, 
Eq. (6.11) of Ref. [10], in the notation of Eq. (4) is 
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This result has a similar frequency dependence as Eq. (8); 
however, its amplitude is much smaller at Ti n 3> T q , i.e., 
the mechanism of Refs. [8-10] appears to be irrelevant. 
Physically, the effect of the ac field on the distribution 
function is dominant because it is accumulated during 
a diffusive process of duration r- m , whereas Refs. [8-10] 
consider only one scattering event. 
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FIG. 1. Photoresistivity (normalized to the dark Drude 
value) for overlapping Landau levels vs lu c /uj at fixed 
ujTq = 2tt. The curves correspond to different levels of mi- 
crowave power 7^ti ^ = {0.24, 0.8, 2.4}. Nonlinear I — V char- 
acteristics at the marked minima are shown in Fig. 2. 

With increasing microwave power, the photoconduc- 
tivity saturates at the value 

Cnh „ ,-9 TtLO TTU) . , ITU) 

= 1 - 8S 2 — cot — , V u sin 2 — > 1. (9) 

Cdc W c LU C U! c 

Note that although the correction is proportional to 
S 2 <C 1, the factor 8n(ui /w c ) cot (ttuj /ui c ) is large in the 
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vicinity of the cyclotron resonance harmonics u> = kuj c 
(k = 1,2,...), and allows the photo- induced correc- 
tion to exceed in magnitude the dark conductivity a^c- 
In particular, cr p h around minima becomes negative at 
Vu > V* > 0, with the threshold value given according 

to Eq. (7) by V* = (V^sin^-sin 2 ^) \ The 
evolution of a B dependence of the photoresistivity p p h 
with increasing microwave power = V u (u c = 0) is 
illustrated in Fig. 1. 

Let us now fix u/ui c such that P* > 0, and consider 
the dependence of <7 p h on the dc field £a c at V u > V*. 
As follows from Eq. (7) , in the limit of large £dc the con- 
ductivity is close to the Drude value and thus positive, 
CTph = (1 — 6<5 2 )er£, > 0. Therefore, cr p h changes sign at 
a certain value £^ c of the dc field, which is determined 
by the condition Q dc = (V u - V*)/V*, see Fig. 2. The 
negative-conductivity state at £d c < £^ c 1S uns table with 
respect to the formation of domains with a spontaneous 
electric field of the magnitude ££ c [7] . 

Using Eqs. (4), we obtain 
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with £* being the threshold value of the ac field at which 
the zero-resistance state develops. Equation (10) relates 
the electric field formed in the domain (measurable by 
local probe [6]) with the excess power of microwave ra- 
diation. It is worth noticing that this relation does not 
include the rate of the inelastic processes. 




FIG. 2. Current-voltage characteristics [ dimensionless 
current j x — (cph/crj^fdc vs dimensionless field Sdc = Q^ 2 ] 
at the points of minima marked by the circles in Fig. 1. The 
arrows show the dc field ££ c in spontaneously formed domains. 

We now turn to the regime of strong B, Ld c T q /ir ^> 1, 
where the LLs get separated. The DOS is then given 
(within the self-consistent Born approximation) by a se- 
quence of semicircles of width 2r = 2(2w c /7rT q ) 1 / 2 : 
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We use Eqs. (1) and (2) to evaluate the OPC at Qd c —> 
to first order in V u and estimate the correction of the 
second order. We obtain 
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The photoresitivity for the case of separated LLs, 
Eq. (12), is shown in Fig. 3 for several values V w of the 
microwave power. Notice that a correction to Eq. (12) 
of second order in V u is still small even at V u > V* = 
T 2 /ujuj c , since lo c P*/T = T/lu <C 1. This means that it 
suffices to keep the linear-in-'Pw term only even for the 
microwave power at which the linear-response resistance 
becomes negative. 




FIG. 3. Photoresistivity (normalized to the dark Drude 
value) for separated Landau levels vs uj c /lj at fixed o;r q = 167T. 
The curves correspond to different levels of microwave power 
= {0.004, 0.02, 0.04}. 

As in the case of overlapping LLs, a negative value 
of the linear-response conductivity signals an instabil- 
ity leading to the formation of domains with the field 
£j c at which cr p h(£dc) = 0. It turns out, however, that 
for separated LLs the kinetic equation in the form of 
Eq. (2) yields zero (rather than expected positive) con- 
ductivity in the limit of strong £dc- This happens be- 
cause elastic impurity scattering between LLs, inclined 
in a strong dc field, is not included in Eq. (2). The inter- 
LL transitions become efficient in dc fields as strong as 
£^ c ~ {Ttr/T c .) 1 / 2 bj 2 /evF [10], which actually gives the 
strength of the field in domains. 

Finally, we calculate the inelastic relaxation time r; n . 
Of particular importance is its T dependence which in 
turn determines that of cr p h. At not too high T, the dom- 
inant mechanism of inelastic scattering is due to electron- 
electron (e-e) collisions. It is worth emphasizing that the 



3 



e-e scattering does not yield relaxation of the total en- 
ergy of the 2DEG and as such cannot establish a steady- 
state dc photoconductivity. That is to say the smearing 
of /o(£) in Eq. (5), which is a measure of the degree of 
heating, is governed by electron-phonon scattering. How- 
ever, the e-e scattering at T » iv c does lead to relaxation 
of the oscillatory term / osc [Eq. (6)] and thus determines 
the T behavior of the oscillatory contribution to cr p h. 

Quantitatively, the effect of electron-electron interac- 
tion is taken into account by replacing the right-hand 
side of Eq. (2) by — St ee {/}, where the collision integral 
St ee {/} is given by 
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and f h (e) = 1 - /(e), e+ = e + E, e'_ = s' - E. The 
function A(E) describes the dependence of the matrix el- 
ement of the screened Coulomb interaction on the trans- 
ferred energy E, 
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where £f is the Fermi energy. Thus A(E) differs from 
the corresponding dependence for a clean 2DEG at zero 
B only by a change in the argument of the logarithm (a 
more detailed discussion will be given elsewhere). 

We linearize the collision integral and solve Eq. (2). 
For overlapping LLs, we put v = 1 in accord with the 
accuracy of Eq. (7) . Then only out-scattering processes 
contribute to the relaxation of the oscillatory part of the 
distribution function (6); the result is obtained by replac- 
ing T in -> T ee (e,T) in Eq. (2) with [12] 
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We turn now to the case of separated LLs. In this 
case, due to oscillation of v, even the linearized collision 
integral gives rise to a non-trivial integral operator. Ana- 
lytical solution of the kinetic equation with this collision 
operator does not seem feasible. However, up to a factor 
of order unity, we can replace the exact collision integral 
with the relaxation-time approximation, thus returning 
to Eq. (2) with 
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One sees that in both cases of overlapping and separated 
LLs the inelastic relaxation rate is proportional to T 2 , 
so that the OPC <7 p h — o"dc m the linear-in-T^ regime 
[Eqs. (8), (12)] scales as T~ 2 . 

Our results are in overall agreement with the exper- 
imental findings [2,3]. The observed T dependence of 
the photoresistivity at maxima compares well with the 



predicted T 2 behavior. Typical parameters u>/2ir ~ 
50 - 100 GHz, r q ~ 10 ps yield lotJ2tt -0.5-1 (over- 
lapping LLs), and the experimental data indeed closely 
resemble Fig. 1. For T — 1 K and ep — 100 K we find 
Tj^ 1 — 10 mK, much less than t^ 1 — IK, as assumed in 
our theory. Finally, for the microwave power — 1 mW 
and the sample area — 1 cm 2 , we estimate the dimension- 
less power <~ 0.005 — 0.1, which agrees with charac- 
teristic values for separated LLs (Fig. 3) but is noticeably 
less than the prediction for overlapping LLs (Fig. 1). The 
reason for this discrepancy remains to be clarified. 

To summarize, we have presented a theory of magne- 
tooscillations in the photoconductivity of a 2DEG. The 
parametrically largest contribution to the effect is gov- 
erned by the microwave-induced change in the distribu- 
tion function. We have analyzed the nonlinearity with 
respect to both the microwave and dc fields. The result 
takes an especially simple form in the regime of overlap- 
ping LLs, Eq. (7). We have shown that the magnitude 
of the effect governed by the inelastic relaxation time in- 
creases as T -2 with lowering temperature. 
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